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Abstract 

Theoretical studies of spin dependent transport in Fe/Cr/MgO/Fe tunnel junctions with non- 
collinear alignment of magnetizations of metallic layers comprising these MTJs are presented. 
Calculations are performed with use of non-equilibrium Green function technique in the 
framework of the Keldysh formalism. The Green and wave functions in the barrier region under 
applied voltage are treated using Wentzel-Kramers-Brillouin (WKB) approximation. Electronic 
band structure of ferromagnetic electrodes is modeled within a two-band model with majority and 
minority states being s-like and J-like electrons, respectively. Furthermore, interfacial s-d 
hybridization is taken into account and calculated using perturbation corrections for the wave and 
Green functions. It is shown that in the presence of Cr layer at the Fe|MgO interface, the 
contribution from s-d hybridization to the total current is much stronger in the anti-parallel 
magnetizations configuration compared to the parallel one leading to decrease of tunnel 
magnetoresistance (TMR) values in agreement with earlier reports. 



1. Introduction 

Magnetic tunnel junctions (MTJ) have been objects of great interest for both scientists and 
engineers because of their high sensitivity to magnetic fields, which makes them good 
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candidates for hard drive read heads and magnetic random access memories ' (MRAM). The 
most attractive among MTJs are their crystalline counterparts due to prediction 3 ' 4 and 
observation ' ' of large tunnel magnetoresistance (TMR) values due to spin filtering of Bloch 
states with certain symmetries in Fe|MgO(001) junctions. It was demonstrated using ab-initio 
calculations 3 ' 4 ' 8 ' 9 that the "s-like" band with A 1 symmetry (existing at the Fermi energy only for 
"up" spins in bcc-Fe (Co and CoFe) along (001) direction) has the smallest decay rate MgO due 
to evanescent state of the same character in this insulator leading to the high conductance in 
parallel (P) magnetizations configuration of adjacent ferromagnetic electrodes. At the same time, 
the conductance in anti-parallel (AP) magnetizations configuration yields very small conductance 
values since there is no "s-like" band for "down" spins in aforementioned bcc ferromagnets but 
"J-like" bands with A 2 , A 2 < and A 5 symmetries. 

Recent exciting experiments 10 demonstrated that additional Cr layer of thickness a inserted 
at the interface between Fe and MgO layers can be used as a tunnel "barrier" for s-like electronic 
states with A 1 symmetry providing straightforward evidence of Bloch states symmetry -based 
spin filtering. In certain sense, bcc-Fe (Co and CoFe) in MgO-based crystalline MTJs should 
demonstrate "half-metallic"- like behavior leading to huge TMR values. However, in 
aforementioned experiment 10 with Fe|Cr|MgO|Fe the TMR strongly decreased as a function of Cr 
thickness and the highest experimental value reported so far in CoFe|MgO|CoFe structures is 
600% at room temperature 11 . One of the possible reasons preventing infinite TMR values can be 
appearance of the partial layer of FeO at the interface which represents an additional barrier for s- 
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electron ' . Another reason could be a symmetry break at the Fe|MgO and/or Cr|MgO interface 
which in turn leads to the s-d hybridization with transition to the conducting state in the anti- 
parallel configuration of magnetizations in Fe electrodes. 



In this work we will theoretically investigate the influence of s-d hybridization on TMR in 
Fe/Cr/MgO/Fe and Fe/MgO/Fe structures with noncollinear alignment of magnetizations of 
ferromagnetic layers comprising these MTJs. 



2. Model 



For sake of simplicity and in a view of arguments stated above, we will describe electronic 
structure of ferromagnetic electrodes within two bands model, i.e. with majority s-like electrons 
and minority J-like holes. It follows from this model that s-electron as well as J-hole bands have 
zero conductivity in the case of the anti-parallel configuration of magnetizations in the outer 
ferromagnetic layers, so the TMR in this structure has to reach infinite value. To calculate the 
conductivity in this non-collinear structure we will use the Keldysh technique for non-equilibrium 
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Green functions with the perturbation corrections coming from s-d hybridization. WKB 
approximation will be used for wave and Green functions in the trapezoidal barrier region under 
applied voltage 14 . To improve further the model we consider different effective mass for electron 
and hole bands in each layer. The details of this approach modified for tunnel transport could be 
found in Ref. [14]. To take into account s-d hybridization one need to find corresponding 
corrections to the wave functions using the following matrix equation: 
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where Vwd)( z )' ^s(d)( z ) are zero 0V ^ QV on hybridization and full wave functions for s, d carriers 



metal/barrier interfacial coordinate (the z axis is perpendicular to the interfaces), and yf a - are 
the parameters of hybridization on the z-th interface. It follows from (1) that the most important is 
assumed to be the hybridization at the metal/insulator interface. In general, there are also terms 
proportional to the off-diagonal on spin parameters of hybridization due to the spin-orbit coupling 
which is neglected here so it's reasonable to assume thaty/ 7 '~ t7 «yf a . All wave and Green 
functions depend on energy E of tunneling carriers and on in-plane momentum k but for 
simplicity we will omit these dependences in expressions below. To find the current we have to 
build the non-equilibrium Green function: 



where indexes /, r mean that that they are the functions for left-right and right-left moving carriers 
correspondently, f(£"), f(£" — eV) - are the Fermi distributions in left, right electrodes, and K is 
the momentum along the interface. Then the current is calculated using the formula: 



and for "up", "down" spins, G°° d ^(z,z{) - retarded Green functions for s, d carriers, z t - is the 



G"+ (z, z') = f(EWi a * (z'Wi a GO + KE - eVWr° * &Wr a GO, 
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The zero order spin "up" s-current and spin "down" d-current are given by the integrals: 
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with corresponding integrands defined as: 
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}LM) = qiqi Re kl j^jT^ ^ [C 1 - cos0 ) Re/C 4 m(\K(q 2 m 2 - ik 5 m b )<pl 

- K~ 1 (q 2 m 2 + ik 5 m b )(pl\ d + (1 + cos 0) Rek$ m 2 \K(q 2 m l — ik^m^cpl 
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where 

K = expQ q(z) dz = exp [— ^ {(gg + k 2 + 0.288 £ ) 3 / 2 - (q 2 + k 2 - 0.288(7 - 
^)) 3/2 }]; ^00 = 2 -^(u-E e -eV^-) + k 2 ; , 1(2) = ,(z 2(3) ) 

In expressions above z 2 , z 3 represent coordinates of the Cr/MgO and MgO/Fe interfaces, 
respectively; Fis voltage, #is the angle between magnetizations of Fe layers, C = O.lSSSmjf 1 x 
10 A/cm ,m 1 ,m 2 - are normalized effective mass for s-electron with "up", "down" spin in Fe 
electrodes, m 3 , m b are the same in Cr and the barrier layers, and 0.288 is factor converging the 
units of energies from eV into A -2 . Furthermore, 

den = (1 + cos 6)[E(q 2 m 1 — ik A m b )cp 2 — £ ,_1 (<7 2 m 1 + ik 4 m b )(pl][E(q 2 m 2 — ik s m b )(p 2 
-E~ l (q 2 m 2 + ik 5 m b )<pi\ + (1 - cos 6)[E(q 2 m 2 - ik s m b )<p\ 

- E~ l (q 2 m 1 + ik 5 m b )(pl][E(q 2 m 1 - ik 4 m b )(p^ - E~ l (q 2 m 1 + ik 4 m b )(pi] 

<pl il) = e ik3a (m 1(2) k 3 - fci (2 )m 3 )(qf 1 m3 + ik 3 m b ) 

+ e- lk 3 a (m 1{2) k 3 + fei(2)m 3 )(<7 1 m3 - ik 3 m b ); 

pt(i) _ e ik 3 a (m 1(2) k 3 - fci( 2 )m 3 )((jf 1 m3 - ik 3 m b ) 

+ e- ik i a (m 1(2) k 3 + /c 1(2) m 3 )(q 1 m 3 + ik 3 m b ); 

where /q, /c 4 (/c 2 , /C5) indicate z-components of s-electron momentum with energy E for "up", 
"down" spins in the left (right) electrodes, and k 3 is the same Cr layer: 



kl { 2) = 7 ^ l iE e -U s ±J sd )-K l ] 

2.m 3 r n 

k! = -^(E e -U^ r )-K 2 . 



Here t/ s - is the bottom of s-band in Fe without exchange splitting J sd , Ug- the same for Cr layer. 
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Fig. 1. Angular dependence of the current in Fe|MgO|Fe tunnel junction with no (jo) and with s-d 
hybridization taken into account. Inset shows a zoom around 9=7t. 

Note that in our model the bottom of electron s-band is situated above the Fermi level giving rise 

to purely imaginary value for k 3 As it follows from Error! Reference source not found, and Error! 

Reference source not found, for purely imaginary k 5 and k% (which is the case in our model) and 

= 7r, both s- and J-currents vanishes giving TMR -> oo. 

In the presence of s-d hybridization additional terms in wave functions appear, for 
example: 

A^ T(T) (z < Zl ) = Griz^yrGW^zdYlWsKzd 

Taking into account additional terms of the second order on s-d hybridization, one gets four new 
currents: two additional ones for s-electron currents due to s-d hybridization at Cr|MgO and 
MgO|Fe interfaces and similarly the two for J-holes: Ajf + AJf T , Ajf + Ajf, Aj$ l + Aj$ l 
and A/3 1 + A/f l . For example, the corrections from the Cr|MgO interface are written as Ajf K + 



A/|k given by the following expressions (we write down here the integrands, and the limits of 
integration were chosen from the condition of reality z-momentum for s- and J-carriers 
respectively): 

T 

+ e ~ ik * a f m d k d + k d m d \}1 . T = ^HL . (6) 
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Finally, the total current is given by the sum of all calculated zero order and additional integrated 
currents: 

h = £ + j& d + A;1 T + A/f + Ajf + A/f T + A;f + A/f + A/ 3 dJL + A/ 3 dJL (7) 
5. Results and discussion 

In order to demonstrate the impact of interfacial s-d hybridization using analytical results of the 
previous section, we calculated the angular dependence of both zero-order and additional currents 
using Eqs (4) and (7), respectively. The results of this calculation for Fe|MgO|Fe and 
Fe|Cr|MgO|Fe structures are shown in Figs. 1 and 2. In Fig. 1 we show the dependence of the 
tunnel current across the standard Fe|MgO|Fe (Cr thickness a=0) junction as a function of angle 
between two Fe magnetizations. As one can see the tunnel current J (filled circles) in the absence 
of s- d hybridization has usual angular dependence and vanishes for antiparallel (AP) alignment of 
magnetizations (Q=n) according to the half-metallic like picture according to our model with 
infinite TMR. At the same time, the total current J-z (open triangles) which takes into account the 
second order contributions due to s-d hybridization is not zero for AP configuration (see inset in 
Fig. 1). This happens because propagating majority (minority) s (d) electrons are converted at the 
interface into d (s) electrons and continues to propagate freely in the layer with opposite direction 
of the magnetization. The same dependences for the junction with finite Cr layer thickness are 
presented in Fig. 2. We notice a strong decrease in the parallel current value since the Cr 
represents an additional barrier for ^-electrons (which is part of Ai Bloch state). In addition, the 
role of s-d hybridization is further enhanced leading to significant drop of TMR ratio. 

We chose the model simulating the half metal nature for both s and J-electrons in the 
absence of hybridization to show more clearly the crucial influence of this hybridization on 
TMR. We have demonstrated that due to the s-d mixing the additional currents comparing with 
the incident s-wp current appear. Due to these processes inside the barrier f d current even 
dominate the /J one because the ip] states easily penetrate it. For the same reason in the region 3 
in the vicinity of the interface the current almost equal to the current because s-electron on the 
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Fig. 2. The same dependencies as in Fig. 1 for Fe|Cr|MgO|Fe tunnel junction with finite Cr layer 

thickness 



interface transform in the d one, which according to the chosen half-metallic model easily 
propagates in the ferromagnetic. 
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